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flj ' Decoherence of a solid state based qubit can be caused by coupling 

to microscopic degrees of freedom in the solid. We lay out a simple 
theory and use it to estimate decoherence for a recently proposed 
superconducting persistent current design. All considered sources of 
decoherence are found to be quite weak, leading to a high quality 
factor for this qubit. 



Q ■ I. INTRODUCTION 

The power of quantum logic ||l| depends on the degree of coherence of the qubit 
C^ . dynamics ^^- The so-called "quality factor" of the qubit, the number of quantum 

^ ' operations performed during the qubit coherence time, should be at least 10^ for 

^^ \ the quantum computer to allow for quantum error correction [y. Decoherence 

is an especially vital issue in solid state qubit designs, due to many kinds of low 
energy excitations in the solid state environment that may couple to qubit states 
and cause dephasing. 

Q^ ' In this article we discuss and estimate some of the main sources of decoher- 

0^ ■ ence in the superconducting persistent current qubit proposed recently H. The 

approach will be presented in a way making it easy to generalize it to other sys- 

2 ! terns. We emphasize those decoherence mechanisms that illustrate this approach, 

and briefly summarize the results of other mechanisms. 

The circuit B consists of three small Josephson junctions which are connected 
in series, forming a loop, as shown in Fig. n . The charging energy of the qubits 
Q ' Ec = e^/2Ci^2 is ~ 100 times smaller than the Josephson energy Ej = ?i/o/2e, 

^ . where Jq is the qubit Josephson critical current. The junctions discussed in H are 

200 nm by 400 nm, and Ej w 200 GHz. 
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FIG. 1. Schematic qubit design [Q| consisting of three Josephson junctions connected 
as shown. Josephson energy of one of the junctions (number 3 in the figure) is adjustable 
by varying the flux in the SQUID loop. The impedances Z^ model electroma gnetic 
environment coupled to the qubit via gate capacitances Cg(i.2). Shunt resistors model 
quasiparticle subgap resistance effect. 



Qubit is realized by two lowest energy states of the system corresponding to 
opposite circulating currents in the loop. The energy splitting of these states 
Eq « 10 GHz is controlled by the external magnetic field flux /, the barrier height 
is ~ 35 GHz and the tunneling amplitude between the two states is i « 1 GHz. 
The Hamiltonian derived in [0 for the two lowest energy levels of the qubit has 
the form 



Tio 



-£0/2 t{qi,q2) 
t*{qi,q2) £0/2 



(1) 



where t{qi,q2) is a periodic function of gate charges qi,2- In the tight binding 
approximation H), t{qi,q2) ~ ti +t2e^"'^^^'^ +^26*^''^^^, where ti is the amplitude 
of tunneling between the nearest energy minima and ^2 is the tunneling between 
the next nearest neighbor minima in the model pj . Both ti and ^2 depend on the 
energy barrier height and width exponentially. With the parameters of our qubit 
design, ^2/^1 < 10"'^, the effect of fluctuations of qi^2 should be small. 

Below we consider a number of decoherence effects which seem to be most rele- 
vant for the design [H, trying to keep the approach general enough, so that it can 
be applied to other designs. 



II. BASIC APPROACH 

We start with a Hamiltonian of a qubit coupled to environmental degrees of 
freedom in the solid: Htotai = Wq((t) +Hbath({fa}), where TLq = Hq + TLcoupUng'- 

nQ = - (A(t) + rfit)) • a , rf^Yl ^""^"" ' (2) 

a 

where <? = {(Tx, <Jy,az) is the vector of Pauli matrices acting on the qubit states, 
the vector A represents external control, and ff is noise due to coupling to the bath 
variables S,a. In (0), A^ = -Eq, A^ - iAy = t{qi,q2). 

The degrees of freedom that may decohere qubit dynamics are: 

• charge fluctuations in the gates coupling qubit states to other states; 

• quasiparticles in the superconductor giving rise to subgap resistance; 

• nuclear spins in the solid creating fluctuating magnetic fields; 

• electromagnetic radiation causing damping of Rabi oscillations; 

• coupling between qubits affecting operation of an individual qubit. 

In all cases except the last one, the qubit is coupled to a macroscopic number 
of degrees of freedom N ^ 1 with about the same strength Aa to each. In 
such a situation, the qubit dccohcrcncc rate is much larger than the characteristic 
individual coupling frequency Aa/h. This means that dephasing happens on a 
shorter time scale than it would have taken to create an entangled state of the 
qubit and one particular element of the bath. In other words, on the decoherence 
time scale each element of the bath remains in its initial state with probability 
1 — 0{l/N), and it is only due to a large number of relevant degrees of freedom 
A^ that the state of the qubit is significantly affected on this time scale. 

This observation makes the analysis quite simple, especially because the condi- 
tion A^ ^ 1 allows one to replace generally noncommuting quantum variables ^q, (i) 
by bosonic fields r]x,y,z(t) fluctuating in time. (Because at large N the commuta- 
tors [?7i(i), '7j(i)] are well approximated by c— numbers.) As a result, the problem 
becomes equivalent to that of longitudinal and transverse spin relaxation times Ti 
and T2 in NMR, corresponding to the noise rji (t) either flipping the qubit spin, or 
contributing a random phase to the qubit states evolution, respectively. Thus we 
can use the standard Debye-Bloch theory of relaxation in two-level systems. 

To adapt this theory to our problem, we assume, without loss of generality, 
that A(i) II z and is constant as a function of time. Then one can eliminate the 
term ^ A • ct by going to the frame rotating around the z— axis with the Larmor 
frequency A = |A|. In the rotating frame the Hamiltonian (|^) becomes: 



^Q = ^ (^11 (0^. + e-^^'v±{ty+ + e'^'vlity^) , (3) 

where 7711 (t) and 'r]±{t) correspond to components of vector r/(t) in (|2|) parallel and 

perpendicular to A, respectively. 

The time evolution due to noise f](t) is given by the evolution operator 

Texp I —i J 7iQ{t')dt' ] written in the rotating Larmor basis. However, for a simple 

estimate below we ignore noncommutativity of different parts of the Hamiltonian 
(|^), and consider a c-number phase factor instead of an operator exponent. 
Then the decoherence can be characterized using the function 

i?(t)-max [(<^|(t)), {\Mt)\^)] , (4) 

where (...) stands for ensemble average, and 

t t 

<P\\{t) ^ j V\\{t')dt' , Mt)=Je-'^''v±it')dt' (5) 



The function R{t) grows with time, and one can take as a measure of decoherence 
the time r for which R{t) ~ 1. There are several assumptions implicit in this 



1) 



which is 
the zero- 



criterion. First, we ignore noncommutativity of different terms in (3 
legitimate at short times, when R(t) ^ 1. Second, we include in (I 
point fluctuations of ^a (t) . The issue of decoherence due to zero-point motion in 
some cases can be subtle. However, since including the zero-point fluctuations in 
R{t) can only overestimate the rate of loosing coherence, one expects the criterion 
R{t) ~ 1 to still give a good lower bound on decoherence time. 

Finally, we note that (|j) contains statistical average over an ensemble of bath 
realizations. Hence care needs to be taken in the interpretation of r when the bath 
is "frozen" into a particular configuration so that the ensemble averaging does not 
apply. In this situation one has to distinguish between decohering individual qubit 
dynamics and averaged dynamics of a qubit array. An example of such a situation 
is provided by the problem of coupling to the nuclear spins, a system with long 
relaxation times. ^ 

Since fj =J2 ^a£,a{t), it is the time evolution of £,a{t) defined by TYbath that is 

a 

what eventually leads to decoherence. One can express quantities of interest in 
terms of the noise spectrum of the components of ff: 

rf^ ' ^^^2 (^ll(-^hlM) (6) 

{\Mm = / du;^-^—^{r^^{~u; ~ A),y^(^ + A)) (7) 



In thermal equilibrium, by virtue of the Fluctuation-Dissipation theorem, the noise 
spectrum in the RHS of (^) and (|^) can be expressed in terms of the out-of-phase 
part of an appropriate susceptibility. 



III. ESTIMATES FOR PARTICULAR MECHANISMS 

Here we discuss the above listed decoherence mechanisms and use the expressions 
(P) and (ph to estimate the corresponding decoherence times. We start with the 
effect of charge fluctuations on the gates due to electromagnetic coupling to 
the environment modeled by an external impedance Z^^ (see Fig. nl), taken below 
to be of order of 400 SI, the vacuum impedance. 

The dependence of the qubit Hamiltonian on the gate charges gi^2 is given 
by dJ), where q\^i vary in time in response to the fluctuations of gate voltages, 
<5<Zi,2 ~ CgbVg(x,2)i where the gate capacitance is much smaller than the junction 
capacitance: Cg <^ C\;2. The gate voltage fluctuations are given by the Nyquist 
formula: (5Vg(-w)JVg(w)) = 2Z^nwcoth Tiw/fcT. 

In our design, |t((7i, (72)! ^ eoi and therefore fluctuations of gi^2 generate primar- 
ily transverse noise 77 j_ in @, ?y±(t) ~ {2'!T/?ie)t2CgSVg{t). In this case, accordingto 
(0) , we are interested in the noise spectrum of SVg shifted by the Larmor frequency 
A. Our typical A ~ 10 GHz is much larger than the temperature ksT/h — 1 GHz 
at T = 50 mK, and thus one has a; ~ A ^ kT/h in the Nyquist formula. 

The Nyquist spectrum is very broad compared to Larmor frequency and other 
relevant frequency scales, and thus in (ITb we can just use the ui = A value of the 
noise power. Evaluating / |(1 — e*"*)/a;f dw — 27rt, we obtain 

2i /27r ^ ^ 



R{t) = {\(^^{t)\') ^ - l—t2Cg j AZ^^A (8) 

Rewriting this expression as R{t) = i/r, we estimate the decoherence time as 



-1 ^ ry-l 

2g2 ^=A \2TrCnt2 



2 



A-^:,3^J=aU-77- (9) 



where h/2e'^ ~ 4kf7. In the qubit design e'^/2Cg ~ 100 GHz, and ^2 - 1 MHz 
when ^2/^1 < lO^'^. With these numbers, one has r — 0.1 s. 

The next effect we consider is dephasing due to quasiparticles on supercon- 
ducting islands. At finite temperature, quasiparticles are thermally activated 
above the superconducting gap Aq, and their density is ~ exp(— Ao/fcT). The 
contribution of quasiparticles to the Joscphson junction dynamics can be modeled 
as a shunt resistor, as shown in Fig. |l|. The corresponding subgap resistance is 
inversely proportional to the quasiparticle density, and thus increases exponen- 
tially at small temperatures: i?qp « i?„exp Ag/fcT, where i?„ is the normal state 



resistance of the junction. For Josephson current /q — 0.2 /iA, i?„ k, 1.3 kfi. At 
low temperatures the subgap resistance is quite high, and thus difficult to measure 
g]. For estimates below we take i?qp = 10^^ 57 which is much smaller than what 
follows from the exponential dependence for T = 50 niK. 

The main effect of the subgap resistance in the shunt resistor model is generat- 
ing normal current fluctuations which couple to the phase on the junction. The 
Hamiltonian describing this effect is 

^cTupUng-Ei^'^^'^W' (10) 



where i = 1, 2, 3 labels Josephson junctions. Projecting ( |10[ ) to the two qubit states, 
one obtains the Hamiltonian (g) with //^(t) = I'l^{t)/e, r]x,y = 0. 

The noise spectrum of the quasiparticle current is given by Nyquist formula: 

(/^P(-tj)/'iP(tj)) = 2R-^hujcoth{hij/kT) (11) 

After rotating the basis and transforming the problem to the form (||) we have 
r]j_{t) ~ {ti/eo)v\\{t), where T]\\{t) ~ I'i^^{t)/e since ti < Eq. 

The analysis of {\(f>±{t)\'^) and (|0||(i)P} is similar to that described 
above for charge fluctuations on the gates, and one obtains R±{t) — 

2t{ti/eo)^hA/{e^Rqp), and i?||(i) = 2tkT/{e^Rqp) which gives 



mm |_rj_,r|| 



e _Rqp f So\ e i?qp 



2hA \tij ' 2kT 



(12) 



Taking Rqp — 10 Ml, T — 50 mK, and So/ti — 100, the decoherence times are 
Til = 1 ms and t±_ — 10 ms. 

The decoherence effect of nuclear spins on the qubit is due to their magnetic 
field flux coupling to the qubit inductance. Alternatively, this coupling can be 
viewed as Zeeman energy of nuclear spins in the magnetic field B(r) due to the 
qubit. The two states of the qubit have opposite currents, and produce magnetic 
field of opposite sign. The corresponding term in (0) is 

^coupling = -Oz E lJ^B{r) ■ s{r) (13) 

r—Ti 

where r^ are positions of nuclei, /x is nuclear magnetic moment and s'(ri) are spin 
operators. 

Nuclei are in thermal equilibrium, and their spin fluctuations can be related 
to the longitudinal relaxation time Ti by the Fluctuation-Dissipation theorem. 
Assuming that different spins are uncorrelated, one has 

{s^{r)s-^{r)) = 2kBT — = ^ ' i^V 

LJ i + U) 1 ^ 



where xo = ^/ksT is static spin susceptibility. 

The spectrum (Oh has a very narrow width set by the long relaxation time 
Ti. This width is much less then ksT and A. As a result, only longitudinal 
fluctuations ?7|[ survive in (^) and (0). One has 



m{t)) = I d^ \^^2j E ^^^B\T){s^{r)s^^{r)) . (15) 

Plugging the spectrum ( p^ in ( [l5| ) and integrating, one obtains 



(16) 



2 \ -1/2 



^0 = ( / %n{r)B^{r)d\ 



where n{r) is the nuclei concentration. The ensembled-averaged decoherence 
time that defined by R{t) ~ 1 is then estimated as: 

To for Ti > To , . 

T^/T, for ri<To ^^'> 

In superconducting Al, nuclear spin relaxation time is strongly varying with tem- 
perature: Ti ~ (300/T[K])e'^/'=«'^s. At T = 50 mK, the time Ti is of order 
of minutes, which exceeds all time scales relevant for qubit operation. To es- 
timate To, we use the magneton /i ~ ehjMc^ where M is proton mass, and 
J B^{r)(Pr ~ 10~^^q/w, where w ~ 0.5 /^m is the thickness of Al wires in the 
circuit, and $0 = hc/2e is the flux quantum. The resulting tq ~ 3 x 10^^ s ^ Ti. 

According to (O), one apparently obtains a worryingly short time t = tq. 
However, we note that this result corresponds to ensemble averaging, and one 
should be careful in applying it to an individual qubit. 

The physical picture is that the nuclei spin configuration stays the same over 
times < Ti . At such times the perturbation ff due to spins has essentially no time 
dependence, and so nuclei can be viewed as sources of random static magnetic 
field. The fluxes of this fleld induced on the qubits depend on initial conditions, 
and are uncorrelated for different qubits. Typical value of this flux corresponds to 
the change in Larmor frequency of order of 5 A — Tq ~ 30 MHz. 

To summarize, for an individual qubit the effect of nuclear spins is equivalent 
to a random detuning caused by random change in A. For an ensemble of qubits, 
there will be a distribution of Larmor frequencies of width SA ~ 30 MHz, even 
if all qubits are identical. However, since the qubit phase can be kept coherent 
within a time < Ti , an indirect observation of Rabi oscillations is still possible by 
using the so-called "spin-echo technique." 



A similar theory can be employed to estimate the effect due to magnetic im- 
purities. The main difference is that for impurity spins the relaxation time Ti is 
typically much shorter than for nuclear spins. If Ti becomes comparable to the 
qubit operation time, the ensemble averaged quantities will describe a real dephas- 
ing of an individual qubit, rather than effects of inhomogeneous broadening, like 
for nuclear spins. 



IV. OTHER MECHANISMS 

Some sources of decoherence are not amenable to the basic approach considered 
above, such as radiation losses which we estimate to have t ~ 10'^ s. 

Another such source of decoherence is caused by the magnetic dipole interaction 
between the qubits. This interaction between qubits is described by 

^coupling = Y. ^^^l""^^ ® ""^^ ' ^^'J ~ I ^'^'|3 (18) 

This interaction is strongest for nearest neighbors. For a square lattice of qubits 
with the spacing R — 10 /im, one has the nearest neighbor coupling A ~ 6 kHz. 
The corresponding decoherence time r = A~^ ~ 0.2 ms is relatively short. 

Several alterations of the design can be implemented to reduce the effect of qubit- 
qubit interaction. One can arrange qubits in pairs with opposite sign of circulating 
currents. This will eliminate dipole moment of a pair, and reduce coupling between 
different pairs to a somewhat weaker quadrupole interaction. The same result can 
be achieved by using a superconducting base plane, in which magnetic dipoles will 
be imaged by dipoles of opposite sign, which will partially cancel the qubit-qubit 
coupling. Also, one can detune Larnior frequencies of neig hboring qubits, moving 



them apart by more than A, which will make couplings (18) off-resonant and reduce 
their effect. 

Unwanted coupling between qubits is a common problem in quantum computers. 
Sophisticated decoupling techniques that have been developed for NMR designs 
g], could equally be used here. The idea is to apply a sequence of single bit 
operations that effectively average out the coupling Hamiltonian over time. Such 
methods would also be effective for reducing the coupling to the environment 0]. 
These techniques are fully compatible with quantum computation and could be 
used to lengthen significantly the effective coherence times. 



V. SUMMARY 

Our analysis shows that for the qubit design g the decoherence time is limited 
by qubit-qubit coupling. By using methods discussed above the decoherence time 



can be made at least 1 ms which for /jtabi = 100 MHz gives a quality factor of 10^, 
passing the criterion for quantum error correction. 

In addition to the effects we discussed, some other decoherence sources are 
worth attention, such as low frequency charge fluctuations resulting from electron 
hopping on impurities in the semiconductor and charge configuration switching 
near the gates [||. These effects cause 1// noise in electron transport, and may 
contribute to decoherence at low frequencies. Also, we left out the effect of the 
ac field coupling the two low energy states of the qubit to higher energy states. 
Results of our numerical simulations of the coupling matrix in the qubit H show 
that Rabi oscillations can be observed even in the presence of the ac excitation 
mixing the states (to be published elsewhere). 
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